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The divergence theorem in its usual form applies only to suitably smooth vector fields. For 
vector fields which are merely piecewise smooth, as is natural at a boundary between regions with 
different physical properties, one must patch together the divergence theorem applied separately in 
each region. We give an elegant derivation of the resulting patchwork divergence theorem which is 
independent of the metric signature in either region, and which is thus valid if the signature changes. 
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(N i I. INTRODUCTION 

In previous work JjJ we discussed the failure, in the presence of signature change, of the standard conservation 
laws normally derived from Einstein's equations. In the process, we derived a form of the divergence theorem which 
is applicable when the signature changes, and related it to the usual divergence theorem across a boundary surface, 
i.e. for piecewise smooth vector fields. We here give a combined treatment of these results which emphasizes their 
similarity. 

After first establishing our notation in Section II, we discuss the divergence theorem in the absence of a metric and 
relate it to the more usual formulation when a metric is given. In Section III we apply the divergence theorem to 
piecewise smooth vector fields, obtaining the patchwork divergence theorem. In Section IV we show how the patchwork 
divergence theorem generalizes standard results on boundary surfaces to our previous results on signature change. 



I Finally, in Section V, we discuss some subtle issues related to the choice of differentiable structure in the presence of 



signature change. 



II. THE USUAL DIVERGENCE THEOREM 



The divergence theorem is usually stated in the presence of a (nondegenerate) metric. However, there is an alternate 
formulation which only requires a volume clement, which we now summarize. 

A volume element on an n-dimensional oriented manifold M is a nowhere vanishing ?i-form which is compatible 
with the orientation on M. In local coordinates x l , an orientation is determined by choosing one of ztdx 1 A ... A dx n , 
and a volume element u> can then be obtained by multiplying this by any strictly positive function ujq. Assuming a 
suitable ordering of the coordinates, we thus have 

oj = ujq dx 1 A ... A dx n (1) 

Given a volume element u> and a vector field X on M, we define the divergence of X by 

div(X) lj := £ x u (2) 

Using the standard expression for the Lie derivative £ x uj of u in terms of components with respect to local coordinates 
x l (see e.g. 0), namely 

(£sW) A ., e = X C d c U a b...e + iOcb...e 9 a X C + LU a c...e ObX" + ... + LO a b...c e X C (3) 

this can be written 
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div{X) = ^3^ +dcX c (4) 

UJQ 



Note that if a (nondegenerate) metric g is given, if u> is assumed to be the metric volume element (so uiq — y] det(g)|) , 
and if V denotes the Levi-Civita connection determined by <?, then our definition agrees with the standard one, namely 

div(X) = V a X a (5) 

This follows since the torsion-free property of V allows us to replace partial derivatives by covariant ones in (@), while 
metric-compatibility means that V a (detg) vanishes. We will refer to this as the physical divergence because of the 
important role played by the metric when making physical measurements. 
Given an open region W of M bounded by S, Stokes' theorem says that 



i= da (6) 
Jw 

for any (suitably smooth) (n — l)-form a, where S must have the orientation induced by W. Using the identity 
relating Lie differentiation to exterior differentiation d and the interior product i, namely (see e.g. g|) 

£ x a = d(i x a) + i x (da) (7) 

for any differential form a, and noting that duo = 0, leads to the following preliminary form of the divergence theorem 



div(X) lj — d> i x uj (8) 
w Js 

In order to relate this to the usual divergence theorem, we need to rewrite the RHS in terms of the induced volume 
element on S. 

We will make the customary identification of the tangent and cotangent spaces of S with the corresponding subspaces 
of those of M, so that in particular we have 

T p S C T p M and T p *S C T* p M (p e S) (9) 

We define a 1-form m 6 T*M for p G S to be normal to S if 

m(Y)=0 {VYeTpS) (10) 

We note that if m ^ is normal to S then each Y G T p M such that m(Y) ^ is not tangent to S, and hence is either 
inward or outward pointing. We further define m to be outward pointing if m(Y) > whenever Y is outward pointing, 
and inward pointing if m(Y) < for all such Y . These definitions extend directly from vectors and 1-forms at a point 
to the corresponding tensor fields. To give a simple and common example, if W — {/ < 0} and S = {f = 0}, then df 
is normal to S and outward pointing. In fact, the 1-forms which are normal to S and outward pointing are precisely 
the positive multiples of df. 

Given a 1-form m which is normal to S, we can define an (n — l)-form a on S via 

m A <7 := u> (11) 

The normality of m ensures that there is a unique a on S satisfying this equation. Since u> is compatible with the 
orientation on M, a is compatible with the induced orientation on S precisely when m is outward pointing. In this 
case, a is the volume element on S induced by m and u>, or more simply the induced volume element on S. The 
interior product is a derivation, so that in particular 

i x (m A a) = (i x m) A a ~ m A (i x a) (12) 

Note that i x m — m{X) and that the pullback of m to S is zero. 

Putting this all together, we finally obtain the divergence theorem in the form 



div(X) uo = <k m(X) a (13) 
Js 

where m is (any) outward pointing normal to <S* and a is the induced volume clement defined above. 
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In the presence of a metric and assuming that S is not null, letting m a denote the (components of the) outward 
unit normal 1-form to S leads to the physical divergence theorem, namely 



V a X a lu= <b X a m a a (14) 
w Js 

where a is the metric volume element on S. It is customary to write this as 

f V a X a d n W = I m a X a d n - x S (15) 
Jw Js 

where d n W = lo and d n ^S = a denote the metric volume elements on W and S, respectively, and m a denotes the 
(components of the) vector field which is the metric dual of m a ; we will refer to this formulation as the physical 
divergence theorem. It is important to note that while m a is indeed a unit vector normal to S, it is outward pointing 
only when it (m a ) is spacelike, and is instead inward pointing where it is timelike; m a is of course always outward 
pointing. 

III. THE PATCHWORK DIVERGENCE THEOREM 

Consider now a boundary surface £, which divides W and S into two parts W + & W~ and S + & S~ , and let 
S° = £ n W be the enclosed region of £. [] Suppose an outward-pointing 1-form m is given on S, and further suppose 
that outward-pointing 1-forms m ± are given on the separate boundaries cW* = S f± U S° of the two regions 
which agree with m on S and which are equal but opposite on E, i.e. 

m \s± = rn (16) 
m~| s = -m + | s =: I (17) 

Let u> and a denote as usual the metric volume elements on M and S, and let a ± denote the induced volume elements 
on the boundaries dW^ 1 , so that 

o- ± \ S =a (18) 
ct-| s = -(7+ls =: a (19) 

We give S° the orientation induced by I and u, namely a, which is the orientation it inherits as part of the boundary 
of W~ (and not W+). 

Now consider a vector field X which is piecewise smooth, so that the usual divergence theorem can be applied in 
each region. Adding the two resulting equations gives 

div(X+) uj+ [ div(X-) uj 
w Jw- 

rn+ (X) cr+ + f m-(X)o-- 

m + (X + )a + -l m + (X + )a + +[ m"(X_)cj"+/ m"(X_) o~ 
Js° Js- Js° 

m{X) <t + + / m(X) <r~ - l([X}) a (20) 
s+ Js- Js° 

where the minus sign in the third line is due to the difference in orientation of S° and dW + and where 

[Q] := lim Q+ - lim Q~ (21) 

denotes the discontinuity in Q across E. We can rewrite this as 



1 We emphasize that E is to be viewed as a hypersurface in a given manifold M, so that there are no complications in the 
manifold structure at or near E; this issue is discussed further in Section 
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/ divpT) ui = I m(X) a - \ l([X]) a (22) 
Jw Js Js° 

and this is the patchwork divergence theorem. In components we obtain 

f V a X a u=l X a m a a- f [X a ] l a a (23) 
Jw Js Js° 

In the presence of a metric, it is customary to assume that m and / are unit; note that the convention adopted here 
is that I points from W~ to W + . If we introduce the (unit) vector field n a which is normal to £ and which points 
from W~ to W + , we finally obtain 

/ V a X a d n W = I m a X a d'^S -el n a [X a ] cf^E (24) 
Jw Js Js° 

where cT 1-1 !] = er° denotes the metric volume element on E and where there is an important sign difference depending 
on whether n a is spacelike (e = I) or timelike (e = —I); l a n a = 1 in both cases. 



IV. APPLICATIONS 



A. Boundary Surface 



A boundary surface in general relativity can be represented as a hypersurface in a Lorentzian manifold across which 
the matter model changes. A spacelike hypersurface corresponds to a change in the matter model at a particular 
time, and a timelike hypersurface corresponds to a change at a particular place, while a null hypersurface corresponds 
to a gravitational shock wave. We thus consider a manifold M with a Lorentzian metric g and a given non-null 
hypersurface £ which divides M into 2 regions M^. Let n denote the unit normal vector to £ which points from M~ 
to M + , and let e = g(n, n) = ±1. The extrinsic curvature K of E can be defined by 

2K := £ n g (25) 



so that in components 

2A' Qb = n c d c g ab + g m bd a n c + g am d b n c (26) 



We shall assume here that g is C 1 across M, so that in particular the Darmois junction conditions 

[h]=0= [K] (27) 

on the induced metric h = g + en®n and extrinsic curvature K of E are satisfied. 

In the non-null case, Israel || used Einstein's equations and the Gauss-Codazzi relations between the curvature of 
M and that of £ to relate the stress-energy tensor of the matter to the intrinsic and extrinsic curvatures of E. Clarke 
and Dray Q generalized some of Israel's results to the null case. Corrected versions of some of Israel's result appear 
in (!]], including 

p := G ab n a n b = X - ((K c c ) 2 - K ab K ab - e TZ) (28) 

where TZ denotes the scalar curvature of the intrinsic metric h on E. In the case where n a is timelike, p can be 
interpreted as the energy density. If the Darmois junction conditions are satisfied, the RHS of this equation is 
continuous at E, and therefore so is the energy density, yielding 



2 We reiterate that we are assuming that E is a hypersurface in a given manifold, so that local coordinates exist which span 
S. Equation (g7|) is thus to be interpreted as applying to n-dimensional tensors, obtained by projection into E. In practice, 
however, it is often convenient to consider the pullback of Equation ( p7| ) to E, which contains the same information but which 
only requires local coordinates within E. 
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[p] = [Gab n a n b ] = (29) 

If we now apply the patchwork divergence theorem (|2^) to 

X a := G ab n b (30) 

then the term containing i([-X"]) vanishes, and we are left with the usual statement of the divergence theorem as though 
no boundary surface were present. If we instead consider 

X a := G a h = hG b a (31) 

we obtain the same result. It is important to note that these results hold even though X may well not be C 1 ! 

The results of Dray and Padmanabhan on piecewise Killing vectors Q can be viewed as application of the patchwork 
divergence theorem in this setting. 



B. Signature Change 



Consider now a situation similar to the above, but where the metric is Lorentzian only in M~, and Riemannian in 
M + , so that the metric is now either discontinuous or zero at E. Such signature-changing models were introduced by 
Dray et al. ; similar models have since been used in a cosmological setting P-|l3| . 

We can impose the Darmois junction conditions by means of 1-sided limits to E. In order to apply the patchwork 
divergence theorem we need a volume element on E, and we cannot now use the metric volume element there. One 
possibility is to work in normal coordinates and note that both sides induce the same volume element; as suggested 
in yj]. We discuss this issue further in the next section, and here assume that a suitable choice has been made. 

On cither manifold- with-boundary M^, let n a denote the unit normal vector as described just before and after 
Even though e now changes sign between the two regions, en a is continuous across E, so that the RHS of (|24) 
remains well-defined. However, Israel's results such as (|2^) no longer hold, and one obtains instead jjj 

[G ab n a n b ] =-U (32) 
[G\ n b l a ] = {K c c f - K ab K ab (33) 

The patchwork divergence theorem now contributes a surface term at E, and for X a = G ab n b as above we now 
obtain Q 

/ V a X a d n W = I m a X a d^S - [ {{K c c f - K ab K ab ) cT^S (34) 
iff Js Js° 

If, however, we instead set X a = G a b l b , we now obtain a different conservation law, namely 

V a X a d n W = <£ m a X a d^S + [ K cf^E (35) 
w Js Js° 

These are two of the main results of @]. 



V. DISCUSSION 



In the normal coordinate approach, M are viewed as disjoint manifolds-with-boundary, with metrics of different 
signatures, which are being identified via an isometry of their boundaries E^. To make the result a manifold, a 
differentiable structure must be specified at the identified boundary E. This can naturally be done by requiring 



3 We will assume throughout this section that these normal vectors make sense, and can be obtained as 1-sided limits from 
M ± . This important issue will be discussed in Section |v], where one way of achieving this is described. As discussed there, 
what really matters is whether the scalars in the last term in ( pi] ) have limits, not whether the normal vectors n a do. 

4 In deriving this result it is slightly easier to use in rather than using directly. 
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that normal coordinates, denned separately on either side of E, be admissible coordinates. In these coordinates, a 
signature-changing metric will take the form 



ds 2 = edr 2 + h tj dx l dx 3 (36) 

where r denotes proper time/proper distance away from E, and {x l : i = 1, ...,n — 1} are local coordinates on E. 
Although this metric is discontinuous at E, the metric volume elements 



d"V := V|dct(.g)| dr A dx 1 A ... A dx™" 1 (37) 

obtained separately on E using 1-sided limits from M ± are identical, so that 

tf^E := v/|det( 5 )| dx 1 A ... A dx"" 1 (38) 

can be taken to be the natural volume element on E. Furthermore, in this approach the normal vectors n a can clearly 
be obtained as 1-sided limits to E. This leads to the results of the previous section. 

On the other hand, one can take the continuous metric approach, in which one starts with a manifold M on which a 
continuous, covariant tensor g of rank two is given which is assumed to be a metric on (the interior of) . The only 
way for the signature of g to change is for it to be degenerate at E; we nonetheless refer to g as a signature- changing 
metric. Such a metric can be put in the form 

ds 2 = Ndt 2 + hij dx l dx j (39) 

where the function N is zero (only) on E. The limit of the metric volume element on M to E is now zero, so that 
it is not clear how to interpret ( pd| ) for the induced volume element. Nonetheless, E does still have a metric volume 
element due to the induced (nondegenerate) metric h, and this volume element agrees with d n_1 E as defined using the 
normal coordinate approach. 

In the continuous metric approach, however, there is no unit normal vector which can play the role of m in the 
divergence theorem. If we pick 

m=^/\N\dt (40) 

so that ( pd| ) is satisfied away from E, then m = on E, and the RHS of the divergence theorem is identically zero 
(provided X is suitably smooth; see below). Contrast this with the normal coordinate approach, in which m = dr is 
a basis 1-form and hence nonzero! In order to preserve a divergence theorem in this approach, we are thus led to pick 
an arbitrary, nonzero m, and then to define the volume element on M via (^) rather than using the metric volume 
element. The patchwork divergence theorem in the form (|2^) now holds, and can be used to derive results analogous 
to those in the previous section, although the volume integral over W is nonstandard. 

Which of these approaches is relevant depends of course on what problem is being discussed. In particular, the very 
requirement that the tensors in the integrands of the divergence theorem be suitably smooth is different in the two 
approaches. For instance, if X a is assumed to be smooth at E in the continuous metric approach, then m a X a = 
in that approach. If, on the other hand, X a is smooth in the normal coordinate approach, m a X a can have a nonzero 
limit at E. In a given physical situation with given smoothness of the physical fields, the same physical results will 
be obtained in either approach (provided suitable limits are taken). However, we feel that, due to the fundamental 
role played by unit normal derivatives in initial-value problems, the normal coordinate approach is more likely to 
automatically incorporate appropriate smoothness conditions at E. 

Clarke & Dray H showed that the continuity of the induced metric h at the identified boundary of two manifolds- 
with-boundary leads to a unique differentiable structure for the resulting manifold. While their derivation assumed 
constant signature, this was only used in showing that the resulting n-dimensional metric is continuous. Thus, their 
result naturally extends to the case of signature change, where it results in the differentiable structure defined by the 
normal coordinate approach (with its discontinuous metric). 

Their result can also be applied to the continuous metric approach as follows. Suppose that one is given two 
manifolds-with-boundary with (suitably controlled) degenerate metrics on the boundaries. Switching to the differen- 
tiable structures on each side induced by (1-sided) normal coordinates leads to a unique differentiable structure via 



5 This problem is in addition to the necessity of specifying which index structure for a given tensor is fundamental, as shown 
e.g. by the difference between Equations (|34[) and (Bq). 
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the generalized result of Clarke & Dray. One can now reverse at least one of the changes in differentiable structure, 
obtaining a manifold structure on the glued-together manifold which agrees with at least one of the manifolds-with- 
boundary. If (and only if) the degeneracies on the two original boundaries are compatible in an appropriate sense 
(based essentially on the behavior of the metric when expanded as a power series in proper time/distance from E), 
a unique manifold structure compatible with both original manifolds-with-boundary is thus obtained. In particular, 
one can always construct several inequivalent manifold structures with continuous (degenerate) metrics, but at most 
one of these will agree with both of the original manifold structures on the two sides. 

This issue can be largely ignored by adopting the invariant approach used in Q, which is closely related to the 
elegant approach used by Carfora and Ellis [jlij. In this approach, one starts by identifying the disjoint manifolds- 
with-boundary M as in the normal coordinate approach. With no further assumptions about the manifold structure 
at E, the 1-sided unit normal vectors still make sense, so that the (pullback of the) Darmois junction conditions can 
still be imposed. 

In this approach, and further assuming that the 1-sided limits of n a X a exist, one obtains the patchwork divergence 
theorem in the form 

/ X7 a X a d n W = I m a X a cT^S - I [en a X a ] cT^E (41) 
Jw Js Js° 

While in practice this is equivalent to (pi}), this form of the theorem emphasizes that it is the limits of physical 
quantities to E which determine the physics, not the choice (or even existence) of a manifold structure there. 

The related issue of defining hypersurface distributions in the presence of signature change will be discussed else- 
where 
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